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Abstract 

In this note, we introduce some models of pedestrian traffic and prove existence and 
uniqueness for these models. 

1 Introduction 

In the last decades, the crowds' dynamics has attracted a lot of scientific interest. A first 
reason of this interest is to understand how crowd disasters happened in some panic events, for 
example at the end of football play, during concerts, in the case of fire, or in place of pilgrimage 
(e.g. on Jamarat Bridge in Saudi Arabia, see |13|). Another reason lies in architecture 
of buildings such as subway stations or stadiums, where a lot of pedestrians are crossing. 
The goals here are consequently twofold: in one hand we want to understand the behavior 
^ _ of pedestrians in panic and adapt the regulation of traffic in order to avoid deaths; in the 

QQ ■ other hand, we want to modelize the interaction of several kinds of pedestrians with different 

. objectives and in particular study how the geometry influences the general pattern. 
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In a macroscopic setting, a population is described by its density p which satisfies the 
conservation law 

dtp + B[v{pVit,x,p))=0, p(0,.) = Po, (1) 

where V{t, x, p) is a vector field describing the velocity of the pedestrians depending on the 
time i > 0, the space x € and the density p. According to the choice of V , various behaviors 
can be observed. Several authors already studied pedestrian traffic in two dimensions space 
' {N = 2). Some of these models are local in p, that is to say V depends on the local density 

p{t, x) [SI [21 [T5l [T6l [20l [21] ; other models use not only the local density p[t, x) but the entire 
distribution of p, for example they depend on the convolution product p{t)*r] |11[I22|. Here, in 
the line of preceding papers 13 [H [3 [6] , we present nonlocal macroscopic models for pedestrian 
traffic, we study these models and compare their properties. 

Our first aim is to modelize the behavior of pedestrians in different situations: crowd 
behaves indeed differently in panic or in a normal situation where courtesy rules do apply. We 
also introduce models in the case of a population interacting with an individual, and in the 
case of several populations with different objectives. For instance, we want to include in our 
study the case of two populations crossing in a corridor. 
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Second, we want to study the introduced models and prove existence and uniqueness of 
solutions under various sets of hypotheses. We will use two kinds of arguments: the first one 
comes from Kruzkov theory \n\ [T8] , the second one from the optimal transport theory. We 
want to prove existence and uniqueness of solutions for the various models presented below. 
Let us concentrate on the case of pedestrians in panic: 

dtp + T)iy{pv{p * ri)v{x)) = 0. (2) 

All the proofs of existence and uniqueness in this note are based on the following idea: let us 
fix the nonlocal term and, instead of we study the Cauchy problem 

dtp + T>i^{pv{r *'q)V{x)) = Q , p(0) = po , (3) 

where r is a given function. Then, we introduce the application 

where the space X has to be chosen so that 

(a) X is equipped with a distance d that makes X complete; 

(b) the application ^ is well-defined: the solution p G X exists and is unique (for a fixed r); 

(c) the application ^ is a contraction. 

Once we have fullfilled these conditions, we can prove existence of a solution using a fixed 
point argument. 

Note that, in the modelization of pedestrian traffic, the space dimension has to be equal 
to two, but our results are in fact true for all N £ N. For instance, they can be adapted in 
dimension = 3 to modelize the behavior of fishes or birds. Consequently, we keep here a 
general N, even if we essentially think to the case N = 2. 

This note is organized as follows: in Section [21 we describe some nonlocal models and their 
properties. In Section [3] we study one of these models through Kruzkov theory and in Section 
[3] we study the same model through optimal transport theory. 



2 Pedestrian Traffic Modefization 

2.1 One-Population model 
2.1.1 Pedestrian in panic 

The first model we present corresponds to pedestrians in panic and was studied in 13 |6] , in 
collaboration with R. M. Colombo and M. Herty. A panic phenomenon appears under special 
circumstances in crowded events. In these cases, the people are no longer rational and try, no 
matter how, to reach their target. Let us denote p{t,x) the density of pedestrians at time t 
and position x € M^. We consider the Cauchy problem: 

dtp + Biv {pv{ri* p{t,x))P{x)) =0; p(0, •)=/Jo- (5) 
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Here, u is a real function describing the speed of the pedetrians. This function does not depend 
on the local density p{t, x) but on the averaged density p{t) * r]{x) = J-^ff p{t, x — y) ri{y) dy. 
The vector field, v{x) describes the direction that the pedestrian located in x will follow, 
independently from the distribution of the pedestrians' density. Note that we are working 
here on the all of and not on a subset of M^; thus, we are not working on a restriction to 
a room, for example. However, we can still introduce the presence of walls and obstacles in 
the choice of the vector field P. Let us denote C the space where the pedestrians are 
authorized to walk, e.g. a room. If we choose v{x) in a nice way (for example we can require 
that on the walls, i.e. for all x G 317, t^{x) coincides with the entering normal to 17), then we 
can conclude to the invariance of the room. More precisely, if the initial density has support 
on some closed set Q, C M^, then the solution will have support contained in O for all time. 
This remark allows us to avoid considering any boundaries and to have solutions on all M.^ . 
Using the Kruzkov theory on classical scalar conservation laws, we are able to prove: 

Theorem 2.1 (see [6]). Let pQ G (L^ nL°° nBV)(]R^, M+). Assume v {'i^'^ , 
V G C^^n W2-i)(M^,M^), T] G C^^p w2,oo-)(-]^Af^]^)_ j^/jg^ ^/^g^g g^^_5^g ^ unique weak entropy 

solution p = Stpo G ^''(M+, L-'-(M^,R+)) to ([5|j with initial condition pQ. Furthermore we have 
the estimate 

||p(t)||Loo < ||po|lL-e^*, (6) 
where the constant C depends on v, v and rj. 

For the definition of weak entropy solutions see Section [3l the proof is defered to Section 
13.2.11 Note that in Theorem 12.11 the hypotheses are very strong. Let us denote 'P(M^) 
the set of probability measures on and A^+(]R^) the set of positive measures on M^. 
In collaboration with G. Crippa, using now some tools from optimal transport theory, we 
obtained the better result: 

Theorem 2.2 (see [9j). Let po G X+(M^). Assume v G (L°° nLip)(M,]R), iJ G (L°° n 
Lip)(R^,M^), T] G (L°° n Lip)(M^,M+). Then there exists a unique weak measure solution 
p G L°°(M4., A^'*"(M^)) to (Ejj with initial condition pQ. 

If furthermore po G L-'"(M^,R+) then for all t > 0, the solution p satisfies also p{t) G 
Li(M^,M+). 

For the definition of weak measure solution see Section the proof is defered to Section 

Note that for the model ([5]), there is a priori no uniform L°° bound on the density. Indeed, 
heuristically, considering the case in which the density is maximal, equal to 1, on the trajectory 
of a pedestrian located in x. If the averaged density around x is strictly less than 1 (because, 
for example there is no one behind this pedestrian), then the speed v{p * rj) will be strictly 
positive, which means the pedestrian in x will try to go forward, even though there is a queue 
in front of him. Consequently, we expect the density to become larger than one. 

This behavior is not really unexpected in the case of panic. In fact, in some events the 
density attained up to 10 persons per square meter, which is obviously too much and a cause 
of deaths (see [13j). Consequently, it is quite satisfactory to recover this behavior. One of our 
goal in this context is then to introduce a cost functional allowing to characterize the cases in 
which the density is too high, and to find extrema of this functional. Let us introduce 

Jt[po) = I I fiStpo)dx , 
Jo Jn 
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where Q C M is the room, po is the initial condition and Stpo is the semi-group generated by 
Theorem 12. II We choose the function / € "^^(M, M+) so that it is equal to zero for any density 
p less than a fixed threshold pc and so that it is stictly increasing on [pc, +00 [. Consequently, 
the functional Jt above allows to characterize the solutions with too high density and in 
particular it vanishes if the set {{t,x) € [0,T] x : Stp{x) > pc} has measure zero. We are 
then interested in finding the minima of this cost functional. 

Using the Kruzkov theory we prove the following differentiability result: 



N 



Theorem 2.3 (see [316]). Let po E (W^'^ n W2'i)(M^, tq G (W^'^ n L~) 

and denote p = StPo- Assume v G ("^^ n W2.~)(M,M), i7 g (-^3 n W2'i)(M^, M^), rj G 
('^^ n W^'°°)(M^, M+). Then there exists a unique weak entropy solution r = S^ro to the 
Cauchy problem 



dtr + Div(r * r]) = —Div{pv'{p * r]) v{x)) , r(0) 



(7) 



Furthermore, the semi-group St obtained in Theorem \2.1\ is Gateaux- differentiable, that is 
to say, for all po G (W^'i n W2'°°)(IR^, tq G (W^'^ nL~^^™^ 



lim 



St{po + hro) - Stpo 



h 



0. 



The proof is defered to Section 13.2.21 

It is not possible to obtain the same result by the use of optimal transport theory. Indeed, 
it seems already not possible to find a good definition of Gateaux differentiability on the set 
of probability measures equipped with the Wasserstein distance of order 1. 



2.1.2 Orderly crowd 

In opposite to the previous model, for a model of orderly crowd it is required to have a uniform 
Jj°° bound on the density. In collaboration with R. M. Colombo and M. Garavello we 
studied the equation 



dtp + Div 



pv{p) v{x 



V{p*r]) 



1+ ||V(/9*7?)| 



(8) 



with Po G (L^ n L~ n B V) (M^ ; M) . 

In this model, the speed v depends on the local density p{t, x), which allows to prove some 
uniform bound in L°°. The prefered direction of the pedestrians is still y{x), but they deviate 
from their optimal path trying to avoid entering regions with higher densities. Indeed, {p * rf) 
is an average of the crowd density around x and —'S/{p * r/) is a vector going in the direction 
opposite to the area of maximal averaged density. Due to the nonlinearity of the flow with 
respect to the local p, it is no longer possible to use optimal transport theory. Hence, we use 
Kruzkov theory, to prove existence and uniqueness of solutions. We obtain the theorem: 

Theorem 2.4 (see [4J). Assume that v G '^^([0, 1], M+) satisfies v{l) = 0, that i7 G (^^ n 



W2,lp^l,oo 



L°°nBV 



an 



N 



d thatr]£ Ctf^ nW^'^ nW^'"^ 



1) . Then, for any po G (L PI 
[0, 1]), there exists a unique weak entropy solution p G '^'^(M+, L''"(M^, [0, 1])). 
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The proof of this theorem rehes on Kruzkov theory (see Section [3]). Note that for this 
model the density is uniformly bounded in L°°, contrarily to the panic model, in which the 
L°° norm can grow exponentially in time. 

Some difficulties now appear in proving that the pedestrians remain in the authorized area. 
Let us introduce the following invariance property: 

(P) Let C be region where the pedestrians are allowed to walk. The model ^ is 
invariant with respect to Q if 

Supp(po) cn Supp(/9(t)) C for aU t > . (9) 



To obtain that (P) is satisfied, we have to require the prefered direction to be strongly entering 
the room (see [4i Proposition 3.1 & Appendix A]). 

Some interesting phenomena show up through numerical computations. First, when con- 
sidering a crowd walking along a corridor, we observe the formation of lanes (see Figure [1]). 



Time = 0.000 



Time = 2.529 



Time = 5.043 




Time ^ 7.557 




Time= 10.071 




Time= 15.014 




Figure 1: Solution to 1^ at times t = 0, 2.529, 5.043, 7.557, 10.071, 15.014. First 3 lanes are 
formed, then the middle lane bifurcates forming the fourth lane. Picture from [Ij. 



Furthermore, this phenomenon seems very stable with respect to initial conditions and 
geometry. Indeed considering the room and initial distribution as in Figure [21 adding some 
various obstacles, we still have lanes, at least in large space (see Figures [3]). 

Time= 0.000 Inside = 100.00% 

ill , , I i 



I 



Figure 2: Initial datum and room geometry. Picture from 



A further remarkable property of the model ([8]) is that it captures the following well- 
known, although sometimes counter intuitive phenomenon (Braess paradox): the evacuation 
time through an exit can be reduced by carefully positioning suitable "obstacles" that direct 
the outfiow (see for instance [12] and the references therein). Indeed, looking at Figure 21 
we observe that the time of exist with obstacles is slightly shorter than the one without any 
obstacle. 
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Time= 2.521 Inside = 76.41% 




Time= 5.043 Inside = 48.41% 
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Time = 7.563 Inside = 24.24% 
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Tinie= 2.521 Inside = 75.20% 



Time = 5.043 Inside = 46.85%. 



Time= 7.563 Inside = 21.41% 
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Figure 3: Solution to jH]) with different geometries, computed at time t = 2.521, 5.043 and 
7.563. Picture from [3]. 



Time= 4.438 Inside = 58.61% 



Time= 6.253 Inside = 41.98% 



Time= 11.396 Inside = 0.99% 
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Time = 4.438 Inside = 59.66% 



Time = 6.253 Inside = 42.96% 



Time= 11.396 Inside = 0.00% 




Figure 4: Solution to ([8]) with e = 0.2, at times t = 4.438, 6.253, 11.396. On the first line, no 
obstacle is present. On the second line, 4 columns direct the crowd flow. The evacuation time 
in the latter case is shorter than in the former one. Picture from [4j. 



2.2 Several populations 

A natural wish now is to extend the previous models to the case of several population with 
different objectives. 

2.2.1 Panic 

With several populations, we have to consider several densities and several equations. For 
two populations, extending the idea of the equation ([5]) to the case of several populations, we 
obtain the system 

dtpi + Div (pi v{pi * 7?i + P2 * m) i^iix)) = , 
dtP2 + 'DiY [p2v{pi*r]i + p2*ri2)iy2(.x)) =0. 

Here we consider that the speed v depends on the average of the total density pi+p2- Further- 
more, the difference of goals of the populations 1 and 2 is reflected in the choice of different 
prefered directions z^i and • 

We are able to prove theorems similar to the ones of section 12.1.11 In collaboration with 
R. M. Colombo |6j, we obtained: 
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Theorem 2.5 (see [6]). Let po = {po,i,Po,2) G (L^ n L°° n BV)(M^,]R^). Assume that, for 

i G {1,2}, Vi e p^2,oo)(]^^j^)^ ^ (^2 p^2,l)(]RAf^^7V-)^ ^. ^ (<^2 p ^2,oo) (jgTV^ j^) _ 

T/ien i/iere exists a unique weak entropy solution p = Stpo G '^''(M+, L-'-(R^, M^)) to (0j with 
initial condition pQ. Furthermore we have the estimate 

IIpWIIlo. < IIPo|lL-e^*, (11) 

where the constant C depends on v, u and rj 

Similarly as for one population, it is possible with several populations to prove the Gateaux- 
differentiability thanks to Kruzkov theory (see [6l Theorem 2.2]). 

As in Theorem 12. 5| the hypotheses of Theorem 12.51 are very strong. In collaboration with 
G. Crippa p)J, using tools from optimal transport theory, we obtained the better result: 

Theorem 2.6 (see [9]). Let po € V{R^f. Assume Vi G (L°° nLip)(M,R), z7j G (L°° fl 
Lip)(M^,M^), r]i G (L°° n Lip)(R^,M+). Then there exists a unique weak measure solution 
p G L~(R+,P(M^)2) to with initial condition pQ. 

If furthermore po G L^(M^,M+) then for all t>0, we have p{t) G L1(M^,M+). 

The idea of the proof for this theorem is given in Section SI 

Interaction continuum / individuals. Note that in the framework of Theorem 12. 6| 
we are dealing with measure solutions. This context allows us to describe a coupling be- 
tween a group of density pi and an individual located in p{t). Indeed, let us assume that 
pi G L°°(R_|_, L-'-(R^, M+)) and p2 = 6p(£^ is a Dirac measure. Then we have the following 
ODE/PDE coupling: 

dtpi + div [pi V {pi * r]i{x) + m{x - pit))) i^i{x)^ = 0, 
p{t) = v(^pi* m (pit)) + r?2(0)) Mpit)) . 



2.2.2 Orderly crowd 

We now extend ([8]) to the case of several populations. We consider here not only that pi and 
P2 have two different goals, but also that pi is repelled by p2, and that p2 is repelled by pi. 
Let us denote the parameter of self-interaction and Eg the parameter of interaction with the 
other population. We obtain: 

dfpi + Divl pi t>i(pi)( i/i(x) - £i -go , II = 0, 

yji + \\vpi*mf \/i + \\'^p2*m\\ 

OtP2+Ulvlp2V2{p2)[j^2[x)-£o , ~ ^' / J J ^ 

^Jl+\\Vpi*r]if ^Jl+\\Vp2*m\\ 

(12) 

Remark 2.7. Note that the interaction between the equations is only in the nonlocal term. 

We obtain again existence and uniqueness of weak entropy solutions thanks to Kruzkov 
theory (see |6l Theorem 3.2]). The obtained theorem is similar to Theorem 12.41 so we don't 
rewrite it. 
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A classical situation considered in the engineering literature is that of two groups of people 
moving in opposite directions and crossing each other. The numerical integration (see Fig- 
ure [5]), shows formation of lanes that are not superimposing as described in the engineering 
literature ([El [10| or [Ij ). 

Time= 0.000 , Time = 1.466 , Time = 1.971 , Time = 2.981 



E 8 -8 -6 



6 8 -8 -6 



Time = 0.000 



Time = 1 .466 





Figure 5: Numerical integration of p2|) . Above, the population pi moving to the right and, 
below, p2 moving to the left. Note the lanes that are formed. First, due to the self-interaction 
{si = 0.3) within each populations and then due to the crossing between the two populations 
(Eo = 0.7). The latter lanes of different populations do not superimpose. Picture from [6]. 

Another situation developping interesting features is the following: two populations are 
initially uniformly distributed in the same region. At time t =0, the first populations starts 
moving towards an exit (first line of Figure [6|) , on the right while the second moves only to 
let the first one pass (second line of Figure [6]). 



Time = 0.000 



Time = 0.203 



Time= 1.213 



Time = 2.476 




Time = 0.000 



Time = 0.203 





-2 2 



Figure 6: Numerical integration of (|12p when with Si = 0, = 0.3. Above, the population 
pi and, below, p2. Note first the formation of small clusters separating the two populations, 
then lanes and, finally, a sort of fingering. Picture from [6]. 



2.2.3 Interaction group / isolated agent 

Let us finally introduce a model that would take into account the situation in which an isolated 
individual interact with a crowd. We think for example to the situation in which a predator 
is running after a group of preys. We can also think to the case of a leader willing to carry a 
group of followers to a given region. Let p € be the density of the group and p € M'^ be 
the position of an isolated agent {e.g. a leader or a predator). We describe the interaction by 
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the coupling (see [7]): 



dtp + Bivi^pV {t,x,p,p{t)) 
p = ip(t,p, {Ap{t)) {p{t)) 



0, 



{t,x) € M+ X 



(13) 



with initial conditions 



p{0,x) = po{x) 



P(0) = Po ■ 



Using once again Kruzkov theory and tools on the stability of ordinary differential equations, 
in collaboration with R. M. Colombo, we proved existence and uniqueness of solutions (see [3 
Theorem 2.2]). 

For example, we can consider: 

• Followers / Leader (see Figure [7|): here the vector p E is the position of the leader 
and p is the density of the group of followers. The function v{p) describes essentially 
the speed of the followers and is, as usual, a decreasing function, vanishing in p = 1; 
the direction of a follower located in x is given by the vector p{t) — x, directed toward 
the leader. The velocity of the leader is increasing with respect to the averaged density 
p*r], computed in the position of the leader. Indeed, we expect the leader to wait for the 
followers to join him when the density of followers is small around him, and to accelerate 
when the density of followers becomes bigger; this means low speed when the averaged 
density is small and high speed when the averaged density is maximal. The direction of 
the leader is a chosen vector field V'(^)- 

dtp + Div (^pv{p) {p{t) - x)e-IIP-^H) = , 
p=il + p*r]{p{t)))^P{t). 



I Solution at time 0,000 



Soiution at time 0.085 




Soiution at time 0,194 




, Soiution at time 0,303 



Soiution at time 0,364 



Soiution at time 0,400 



Figure 7: Numerical solution of (|13p in a predator/prey interaction. Picture from [7] 



• Sheeps / Dogs: in this example, there are n dogs, located in Pi(t) G M? for i G {1, . . . , n} 
and a group of sheeps of density p{t,x). The dogs want to constrain the sheeps to a 
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given area by running around them and the sheeps are afraid by the dogs. As above, 
the speed of the sheeps is given by the decreasing function v{p) that vanishes in p = 1. 
The direction of a sheep located in a; is a sum of two terms. The first one is their 
prefered direction I'ix) ; the second one is given by the vector Yl^=ii^ — Pi{t))e~^^'P^~^^^ 
representing the repulsive effects of the dogs on the sheeps, and its presence is is less felt 
if it is too far away. Each dog run around the flock tacking the direction perpendicular 
to the direction of maximal averaged density. 

dtp + Bw (^pv{p) (i7(x) + Er=i(3;-pOe""^'""")) =0, 
{p*Vv)^{Pi{t)) 



Pi 



l + \\p*Vri{p,{t))\ 



for all i € {1, . . .n} . 



• Preys / Predator: here a predator, located in p{t) € M'^ is pursuing a group of preys of 
density p{t,x). For example, we can think to an hawk pursuing a swarm of little birds 
or to a shark pursuing a shoal of little fishes. As above, the speed of the preys is given 
by the decreasing function v{p) that vanishes in p = 1. The direction of a prey located 
in X is given by the vector (x — p(t))e~ll^^*^~^ll so that the predator has a repulsive effect 
on the preys, and its presence is is less felt if it is too far away. The acceleration of 
the predator is directed toward the maximal averaged density of preys, as felt from its 
position, so it is directed as V(p * r]){p{t)). 

dtp + DiY(^pv{p) (l+ e-il"-P(*)ll (x-p(t)))^ =0, 
^ = ip* V)r?(p(t)) . 



Note that Theorem 12.61 deals with measure solutions. Consequently, if pi G L-*" and 
P2 = ^p{t)i ^6 recover the above coupling PDE/ODE. However, in Theorem 12. 6| there are 
restricitions on the fiow as we only deal with continuity equations. In particular, with a well- 
chosen nonlinearity with respect to p in the first equation of (jl3p . we are able to provide for 
()13p an a priori uniform bound on the L°° norm, which is not possible for (|10p : in the three 
examples above, if f (1) = 0, then < po ^ 1 implies < p{t) < 1 for all t > 0. 



3 Using Kruzkov theory 

3.1 General Theory on scalar conservation laws 

In this section, we consider classical scalar conservation laws. 

dtu + Dwf{t,x,u) = F{t,x,u), u{0,x) = uo{x) , (14) 

where / is the flow and F is the source. Note that, Div stands for the total divergence 
whereas div stands for the partial divergence. Thus 'Divf(t,x,u{t,x)) = divf{t,x,u{t,x)) + 
dufit, X, u{t, x)) ■ Vu{t, x). 

Let us recall the deflnition of weak entropy solution: 
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Definition 3.1. The function p € L°°([0, T] x M^,M) is called weak entropy solution to 
the Cauchy problem (|14p if it satisfies for all A; E M and any test-function Lp € 'io^(\ — 



cx),r[xM^ 



(15) 



(n - A;) dt^p + {f{t, x, u) - f{t, x, k)) V^ip + x, u) - div/(t, x, k)) ip 

xsign(ti — k) dx dt + \pa{x) — k\ip{<d, x) dx > . 
We also recall the well-known Kruzkov theorem |17( Theorem 4], giving existence and 



1). For all yl,r > 0, we denote = 
Under the conditions f G '^°(0;M^), F G 



uniqueness of weak entropy solutions. 

Theorem 3.2 (Kruzkov). Let uq G (L°° n 
[0, r] X X [-A, A] and n = R+ xR'^ xM 
•rO(ri;M) and 

' / G '^°(J7;M^) , FG'rO(n;M), and for all T, A > : 

f , F have continuous derivatives: duf , duVf , V^/ , duF , \/F , 
for aUT,A>0 , duf G L~(Q^) , 
F-div/ G L°°(fi^), du{F -divf) e'L°°{n^) 



(K) 



of J j^| ) i/iai is right- 



there exists a unique weak entropy solution u G L°°([0, T]; L-*- 
continuous in time. 

Let vq G (L-*- n L°°)(M^;IR). Let u he the solution associated to the initial condition 
uq and V he the solution associated to the initial condition vq. Let M he such that M > 
sup(||«||l^(k+xRJv.r), ||f||L^([o,T]xiRJv.R)). Then, for all t G [0,r], with j 



have 



Ili < e"'*\\uo 



voW-Li 



3.1.1 Estimate on the dependence with respect to fiov^ and source 



Let us now consider a pair of difi^erent flows f,g £ '^^([0, T] x 



(16) 



and different 



sources F,G £ '^^([0,r] x x R;R), let us denote u,v the solutions of 

dtu + Div f{t,x,u) = F{t,x,u) , u{0,-)=uq, (17) 

dtv + Div g{t, x, v) = G{t, x,v), v{0,-) =vo, (18) 

with initial conditions Mo, fo G (L-'-nL°°nBV)(M^, M). We want here to estimate || (n — f)(t)||j^i 
with f—g, F — G, Uo — vq. To do that, we use the doubling variable method due to Kruzkov [17J. 
This strategy was already employed by Lucier |19) and Bouchut & Perthame [3] to study the 
case in which the flows /, g depend only on u and the source are identically zero {F = G = 0). 
A key ingredient in the previous results is that the solution of dtU + Div{f{u)) = with initial 
condition uo G (L°° n n BV)(M^, M) satisfles TY{u{t)) < TV(no). 
Let us remind the deflnition of total variation. 



Definition 3.3. For u G L/^^ 



we denote the total variation of u: 



TV(n) 



sup 



udiv^; ^ £'€l[ 



I^IIl- <i 



The space of function with bounded variation is then deflned as 



BV( 



u G L/-„^;TV(n) < ooj 
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Remark 3.4. If u G (^^ n Wi'i)(M^, M) then TV(n) = |IVu|Ili. When / and F are not 
depending on u, we have 

no G (L~ nBV)(M^,M) ^ Vt > , u{t) G (L°° nBV)^™^ 

and, with 7 = ||i9„F||loo(qM), 

TV('u(t)) < TV(no)e^* . 



Theorem 3.5 (see [19]). Letf,g:R^ be globally Lipschitz, letuo,Vo G {1,^(11^°° 
be initial conditions of 

dtu + Biv f{u) = 0, dtv + DiYg{v) = 0. 

Assume furthermore that vq G BV(]R^;]R). Then for all t >0, 

\\{u - v){t)\\.^i < \\uo -voIIli +CtTY{vo) Lip (/-fir), 

where C > is a given constant. 

We want to generahze this theorem. Let us introduce the set of hypotheses: 

' for all [/, T > , duF G L°°(f^^) 
(FS)i fT r 

Denoting Q'f = [0,T] x x [-y,^], we obtain: 



Theorem 3.6 (see [IS]). Assume that {f,F),{g,G) satisfy (K), that {f — g, F — G) satisfies 
(FS). Let Uo,Vo G (L^ n L°° n BV)(R^;R). Let u and v be the solutions of fl^) and 
associated to {f,F) and {g,G) with initial conditions uq and vq. Assume furthermore that 
TV(n(t)) < 00 for all t > 0. 

Let V = max(||u||Loo, ||'v||lcx)) and k = \\duF\\^^^^vy Then for all t G [0,T]; 



I, J,; 

ft 



\{u-v){t)\\^^<e^'\\uo-vo\\-^i+e''' sup (TV(u(r))) / \\duif - 9)iT)\\^o.^^N ^^w])^'^ 

rG[0,t] ^0 \ I ^ I) 



+ Te^^*"") / ||((F-G)-div(/-5))(r,x,-)| 

Jo JR^ 



L°°{l-V,V])^^^'^ ■ 



The proof of this theorem is based on the Kruzkov doubling variables method and is 
detailed in [18j. 

In some particular cases, we recover known estimates: 

• f{u),g{u),F = G = : k = and 

||(u - w)(t)||j^i < ||no-uo|lLi +*TV(mo)||5„(/ -5)||Loo(nV) • 



f{t,x),F{t,x) : K = 0, u{t,x) = uq{x) + (F — div/)(r, x) dr, and 

\\{u-v){t)\\^,<\\uo-vo\\i^i+ f j |((F-G)-div(/-5))(T,x)|dxdr. 

Jo Jr^ 
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3.1.2 Total variation estimate 

To complete Theorem [321 we need an estimate on the total variation, defined in Definition [3]3j 
To do so we need a few more hypotheses. Indeed, there is no reason why the total variation 
should be bounded for all time. In fact, let us consider the equation dtu + dx cosx = 0, uq = 0. 
The solution of this Cauchy problem is u{t,x) = tsinx whose total variation is at time 
and +00 for any time t > 0. 

Our goal now is to give a general estimate on the total variation. Let us introduce the 
following set of hypotheses 

' for all A,T >0 

V(F-div/)(t,x,-)||Loo([_^,^].KiV) < oo. 
We obtain 

Theorem 3.7. Let us assume {f,F) satisfies (K)-(TV). Denote Wn = /J'^^(cos (9)^ de*, 
M = ||ii||l°o([o,t]xRJV); and 

KQ = {2N + l)||Va„/||Loo(f^M) + ||9«i^||Loo(f^M) . 

Then, the weak entropy solution u of J j^p satisfies u{t) € BV(M^;M) for all t > 0, and 
TV(n(r)) < TVMe'^°^ + || V(F - div/)(t, x, OlLoofr ^z, c/,!*^^ , 

where Ut = supy^^N \u{t,y)\. 

The proof of this theorem is based on a good representation formula for the total variation 
and on the doubling variable method and is detailed in )18j . 

Remark 3.8. In some cases, we recover known estimates. 

• When / depends only on u and F = 0, we have a result similar to the one that was 
already known: TV(n(t)) < TV(no). 

• When / and F do not depend on u, the equation reduces in fact to the ODE dtu = 
(F — div/)(t, x), whose solution writes 

u{t,x) = uo{x) + / (F — div/)(T, x) dr . 
Jo 

Meanwhile, the bound above reduces to 

TY {u{t)) <TY{uo) + NWn [ [ |(F-div/)(r,x)| dr 

Jo JR^ 

which is essentially what we expected, up to the coefficient NWn- 
Proof. Let us admit the following proposition: 
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Proposition 3.9. Let ^ G '^c°°(]R+; M+) such that ||//||li = 1 and fi' < on M!j_. Let 
u € LJ-^^(IR^,M). We define fJ-x{x) = pvA* ■ ^^^'"^ exists Cq > such that VA > 0, 

^W = t[ I \u{x + y) -u{x)\fix{y)dxdy <Co, 



then u G BV( 



and 



CiTV(u) = limX(A) < Co. 

A->0 



with Ci = /jgjv \yi\ti{\\y\\)dy. 
Hence, we introduce: 



T{T, A) = / / / 



R+M{To-t)) 



\u{x + y)- u{x)\fixiy)dx dy dt , 



our goal being to estimate tliis quantity using tlie Kruzkov doubling variable method. 

Doubling variable method. Let us denote u = u{t, x) and v = u(s, y) for (t, x), {s, y) € 

Ml X R^. For all k,l and all test-function (p = (p{t,x,s,y) £ 'tf^ 



have, by definition of weak entropy solution (see Definition 13. ip 

(n - k) dtip + (/(t, X, u) - f{t, X, k)) V ^^p + x, u) - div/(t, x, A;)) tp 

xsign('u — fc) dx dt > 



we 



(19) 



(20) 



and 

/ / {v-l)ds^+{f{s,y,v)- f{s,y,l))Vyip + {F{s,y,v)- d:\Yf\s, y,l))Lp 

xsign(f — /) dy ds > . 

Let us choose k = v{s, y) in (|19p and integrate in (s, y). Similarly, we choose / = u{t, x) in (|20p 
and we integrate in (t, x). Furthermore, we choose (p(t, x, s, y) = ^{t — s,x — y)^{t, x) and we 
summate to get 



/ / / / sign(u - v) 

JR* JM'V Jjj* JjjAT 



(u-v)^ dt*^ + {fit, X, u) - fit, X, v)) ■ (V$) ^ 



+ if is, y, v) - fis, y, u) - fit, x, v) + fit, x, u)) ■ (V^) ^ 



(21) 



+ {Fit, X, u) - Fis, y, v) + div/(s, y, u) - div/(t, x, v)) if 



dx dt dy ds > . 



We choose now /i and v unit approximations and 

^it — s,x — y) 



A^ 



\ \ J ri \ ri 



so that ^ converges to a Dirac in space and time when A, r/ — )■ 0; and we choose <1> such that 
$ — )• l[o,r]x_B(xo,/?+A/(T-t)) when 0, e — > 0. The last step consists in estimating the integral 
()2ip when rj,9,£ — t- 0, keeping the A fixed to make appear J^. 
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Estimate on T . By the doubling variable method, we obtain 

5r-F(r, A) <aT-F(0, A) + C\ (dxT{T, A) + y ^(T, A)^ + C'T{T, + >^ £ , (22) 

where 

A{t) = Mil ||V(F-div/)(t,x,.)||Loof„)dx, Mi=/ \\yM\\y\\)dy, 

Jm ^ ' Jrn 

C = A^||V9„/||loo + \\duF\\^^ . C = ||V9„/||l.. . 

We use the hypothesis uq € BV(M^, M) to get ^9t-F(0, A) < MiTV(no), and we integrate 
with respect to time, obtaining: 

< ^TV(no) + dxHT, A) + '^T{T, ^) + ^ £ Mt)dt , 

where 

C 1 

a{T) = N + — - — ^T^o -oo . 
Next, we choose T so that a < —1 and we integrate on [A, +oo[. We obtain: 

/•(T,A)<^-^TV(no)+ ^ rA{t)dt. 
—a — 1 C C[—a — i) Jq 

Besides, note that, thanks to the shape of ^, we can find > so that 

dxHTA)+^J' = r I I \u{x + y)-u{x)\i^'x{y)dxdydt < ^T{T,2X) . 

^ Jo JRJV J B{xo,R+M{To-t)) 

Hence, we finally obtain 

\dTT{T, A) < ^5r/-(0, A) + ^^^^^^^ {m{Y:N(u^) + ^'^ A(t)dt\^ + A{t)dt . 

By Proposition 13.91 we have then uit) S BV(M^,M).We can now divide ([22]) by A and make 
A — 7> 0. Applying a standard Gronwall-type argument, we obtain the desired inequality. □ 



3.2 Scalar conservation law with a non-local flow 

Let us consider equation ([5]) with initial condition po ^ (L^ n L°° n BV)(M^,M). We 
want to prove Theorem 12.11 and Theorem 12.31 Note that this proof adapt smoothly to 
prove Theorem 12.41 and 12.51 We use here the scheme proposed in Section [H choosing X = 
^°([0,r],L^(M^,]R)) at point (a). Point (b) is then satisfied by Kruzkov theorem (see The- 
orem 13. 2p and point (c) is satisfied by the stability estimate given by Theorem 13.61 & 13.71 

This scheme of proof can be also used for the other models: the proof of Theorem 12.41 for 
an orderly crowd is based on it, as well as Theorem 12.51 giving existence and uniqueness for 
the several populations models (jlOp and (|12p . For more details on these various theorems, 
refer to [6lll]. 

We use also this scheme of proof in order to prove existence and uniqueness of solution to 
equation ()13p . where we need furthermore stability results on the ODE (see [7]). 
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3.2.1 Proof of Theorem [2711 



Let us introduce the space 

X = <r°([0,r],Li(M^;M+)). 

Let r E X and p E X be the solution of ([3]) given by Kruzkov Theorem (see Theorem I3.2p . 
Note that, u = being a solution to ([3]), /Oq > implies p{t) > for all t > thanks to the 
maximal principle on scalar conservation laws |17| Theorem 3]. 

Remark 3.10. For the model ([8]) of orderly crowd, we can use the same king of argument to 
have a uniform bound on the L°° norm. Indeed, if u = 1 is a solution to the equation 

dtp + DiY {pv{p)w{x)) = (23) 

then po < 1 implies p(t) < 1 for all t > 0. This fact allows us in this case to consider the space 
X = '(r°([0,r],Li(R^, [0, 1])). The same fact is also true for each density of the model ([12]) 
of orderly crowd with several populations. Indeed, the interaction between the equations is 
only in the non-local term; so, if we fix the nonlocal term, we have to deal with the same kind 
of scalar conservation law as ()23p . In particular, for each density pi the maximum principle 
gives us 

< po,i < 1 < pi{t) < 1 for aU t > . 

So finally, if there are k populations, we have the following uniform bound on the total density 
density 

k 

0<^Pi{t) <k. 
1=1 

Let us introduce the application 

^■.reX^peX. 
For ri , r2 , we get by Theorem 13.61 and Theorem 13. 7| 

ll^(ri) - ^(^2)||loo([o,t],L1) ^ /(^)ll^i - '•2|Il-{[0,T],L1) ' 

where / is continuous, increasing, /(O) = and / —j't^oo oo. Consequently, we can choose T 
small enough so that f{T) = 1/2 and apply the Banach fixed point theorem for a small time. 
Iterating the process in time, we obtain global in time existence. 



3.2.2 Gateaux derivative of the semi-group 



We prove now that the semi-group obtained in Theorem 12.11 (or 12. 5p is Gateaux-differentiable 
with respect to the initial condition. The following result is obtained by the use of Kruzkov 
theory. 

Note that the Gateaux-differentiability is also possible for the more non-linear models of 
orderly crowd jS]) and (jl2p . Indeed an important tool below is that the equation ([5]) preserves 
the regularity of the initial condition, which is no longer true for ([8]). 



Definition 3.11. We say the application S : L^{R^;R) L^{]1 



is Gateaux differen- 



tiable in li^ , in po E L-*-, in the direction tq E L-*- if there exists a linear continuous application 



DS{po) : 



L-*- such that 

S{po + hro) 



S{po 



h 



DSipo)iro) 
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We expect the Gateaux derivative to be the solution of the hnearized problem ([7]) . 
First, we prove that the linearized problem admits an entropic solution (see 

Theorem 3.12. Assume that v £ (<^4 n W2.~)(M, M), P e (<^3 p ^2,1 p ^2,oo^mAf ^^.iv^ 



r] G (<(f3nw2'°°)(M^,M+). Letpe '^^{[0,Tecc[;W^'°°nW^'^{R^ ,R)), tq G (LinL°°)(M^; M). 
Then the linearized problem ^ with initial condition tq admits a unique entropic solution 
r G '^°([0,re^[;Li(M^;]R)) and we denote Y^^ro = r{t, ■). 

If furthermore tq G W1'1(R^,M), then for all t G [0,Tea;[, r{t) G Wi'i(R^,M). 

The proof is similar to the one of Theorem 12.11 so we omit it. 

Theorem 3.13. Assume that v G n W2'~)(M, M), P g (<^3 p ^2,1 n w2''^)(M^, M^), 
rj G ("^3 n w2'~)(M^,M+). Let po G (W^'"" n W2'i)(M^, M+), tq G (W^-^ n L°^)(M^, M). 

Then, for every t > 0, the local semigroup of the pedestrian traffic problem is L-*- Gateaux 
differentiate in the direction tq and 

DSt{po){ro) = Yf'''ro. 

We prove this theorem following the same sketch of proof as above. 

Proof. Let p,ph the solutions of the problem jS]) with initial condition po, po + hrQ. 

Let r be the solution of the linearized problem ([7]), r(0) = ro and let Zh = p + hr which 
then satisifies 

dtZfi+Di-v (^Zh{v{p * rj) + hv'{p * r]){r * r/))z7(x)) = h'^DiY{rv' {p*r]){r*ri)u{x)) , Zh{0) = po+hr^ . 
Then we use Theorem 13.61 fc [3^1 to compare ph et Zh. We get: 

+ /iC(/3)Te^(«^||r||Loo(wi.i)lkllLoo(Li), 

with F increasing and -F(O) = 0. Then, we choose T small enough so that F{T) = 1/2. 

Besides, applying Theorem 13. 6l to compare p and ph and using Gronwall lemma, we obtain 
that j^Wph ~ pIIl°°{L1) remains bounded when /i — > 0. 

Finally, we can take the limit h ^ and obtain the desired result. □ 



3.2.3 Extrema of a cost functional 

We are now able to characterize the extrema of a given cost functional. Let J be a cost 
functional such that 

JiPo)= f (StPo) i'it,x)dx . 



Proposition 3.14. Lei / G '^^'^(M; M+) and ^ G L~(M+ x M^; M). Assume S : [0,T]x{L^n 
L~)(M^;E) ^ (L^ nL~)(M^;M) is Gateaux differentiable. 
If Po G (L^ nL°°)(M^;M) is solution of the problem: 

Find min J(/9o) such that\StPo is solution of problem (Traffic)} . 

PO 

then, for all ro G (L^ n L°°)(M^;M) 

/'(5tpo)SfroV(t,x)dx = 0. 
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4 Using the Optimal transport theory 



In this section, we want to prove Theorem [221 giving existence and uniqueness of weak measure 
solurtion to ([5]) and (jlOp . In a more general settings, we consider the system: 



dtPi + d\Y{piVi{x,pi *r]i^i,...,pk *rii,k)) = 0, 

Pi(0) = Pi , 



(t, x) gM+ X (24) 

i e{i,...,k}. 



in which the couphng between the equations is only present through the nonlocal term. Fixing 
the nonlocal term, we obtain a system of decoupled continuity equations: 

dtpi + div(pi bi{t, x)) = 0, 
dtPk + div(/3fc hk{t, x)) = , 
where 6i, . . . , 6^ are regular with respect to x. 

4.1 Existence and uniqueness of weak measure solutions 

Let us remind the following definitions of weak measure solution, but also of Lagrangian 
solution. 

Definition 4.1. p G L°°([0, T], A^+(M^)'=) is a measure solution of ([MD if, G "^^i] - 
oo,T] X 



JR^ '- 



dtV + Vi{x, p*T]i) ■ Vip dpl{x) dt + 



/ ip{0,x)dpi{ 



x)=0. 



Definition 4.2. p G L°°([0, T], is a Lagrangian solution with initial condition 

p G 7W+(M^)*-' if there exists an ODE flow X' : [0, T]xR^ ^ M^, solution of 

dX 



it,x)= Vi{X'{t,x),pt*r,\X'{t,x))), 
X'{0,x) 



and such that pi = XI ^p^ where XI : 
any {t,x) G M+ x M^. 



x; 

i>N , mN 



is the map defined as Xl{x) = X^{t,x) for 



Definition 4.3. Let p be measure on Q and T : Q Q' a measurable map. Then T^p is the 
push-forward of p if for any ip G 



/ ip{x) dT^p{x) = j if {T{y)) dp{y) . 
Jn' Jn 



If Tjj d^ (x) = /(x) dx and d/i (y) = dy; and T is a '^^-diffeomorphism. Then we have 
the change of variable formula 

5(x) = /(r(x))|det(vr(x))|. 

Instead of proving directly Theorem 12.21 (or Theorem 12. 6|) . we prove first the existence and 
uniqueness of lagrangian solutions. 
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Theorem 4.4. Let p G A^+(M^,M'=). Let us assume that V G (L°° n Lip)(M^ x M*=,M^x'=) 
and that r] G (L°° n Lip) (M^,M*= Then there exists a unique Lagrangian solution to system 
^24\ ) with initial condition p. 

Remark 4.5. Note that a Lagrangian solution is also a weak measure solution. 

Besides, we recover the results of Theorem 12. H except the continuity in time. Indeed, 
Lagrangian solutions conserve the regularity of the initial condition along the time: if p G 
Li(M^, (M+)'') then p G L°°(E+, L1(M^, (M+)'')) and we have 

If moreover ^ G (L^ n L°°)(M^, (M+)^), then p{t) G L°° for all t > and we have the 
estimate 

IIpWIIl- < llp|lL-e^S 

where C depends on V and ry. 

We conclude the proof of Theorem 12.21 stating the uniqueness of measure solutions, so that 
the Lagrangian solution obtained is also the unique measure solutions. 



Proposition 4.6. Under the same set of hypotheses as Theorem \4-4\ ^he weak measure solu- 
tions for are unique. 

4.2 Proof of Theorem 14.41 

Once again, the proof is based on the scheme described in Section [TJ We consider below 
probability measures instead ; for any i, we could also consider positive measures of fixed 
total mass p* 

M 

(a) Let us first introduce a space 



N\k\ 



X = L°°([0,r],P(M 

We equip this space with the distance d{p, u) = sup^gjOjT] ^t)- The Wasserstein distance 

Wi is defined as follow: 

Definition 4.7. Let p,v & P(M^). Let us denote : M'^xM'' the projection on the first 

coordinate; that is, for any (n, f ) G M'^ x W^, Px(i*, v) = u. In a similar way, Pj^ : R'' x M'^ — ^ M'^ 
is the projection on the second coordinate; that is, for any {u,v) G M*^ x M*^, ¥y{u,v) = v. We 
denote S (p, i/) the set of plans, that is 

H(p, I/) = {7 G V{MJ^ X M^) : P^0 = p and ¥y^^ = . 
The Wasserstein distance of order one between p and v is 

Wi{p,v)= inf / \x-y\d"f{x,y) . 

Let p = {p^, . . . , p^), a = {a^ , . . . , a^) G V{M.^)^ . The Wasserstein distance of order one 
between p and cr, denoted Wi(/0, cr), as 

k 

Wi{p,a)=Y,Wi{p\a'). 
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Let us recall the following duality formula: 

Proposition 4.8 (cf. Villani |23l p. 207]). Let f,g be two probability measures. The Wasser- 
stein distance of order one between f and g satisfies 

Wiif,g)= sup / ip{x) {df{x) - dg{x)) . 

Lip((/p)<l JR'* 

(b) Existence of Lagrangian solutions. 
Let r G L°°{[0,T],V(R^ f). 

Define bi{t,x) = Vi{x,rt * r]i) G L°°([0, T], Wi'~(M^)'=). Let us consider the equation 

dtpi + div {pibi{t, x)) =0. (25) 
Let Pi = Xt^p he the Lagrangian solution of (j25p . Then the application 

: r G L~([0,r],P(IR^)^') ^ /> G L°°([0, T], P(M^)'=) . 

is well-defined. 

(c) Stability estimate. 

Proposition 4.9. Let p,a eViR^) and r,s e '^^{[0,T],'P{M.^)). Lei F G (L~ n Lip)(]R^ x 
M'=,M^), r],ue (L°°nLip)(M^,M). If p and a are Lagrangian solutions of 

dtp + div{pV{x,r*T]))=0, p(0, •)=^, 
dta + div((T V{x, s * rj)) = , a{0, ■) = a . 

We have the estimate: 

Wi{pT,aT) <e^^Wi{p,a)+Te^^C' sup Wi{rt,st), 

te[o,T] 

where C = Lip^,(F)+Lip^(y)Lip(r/)||p||_^+Lip,,(F)Lip(7/)|j^||_^ andC = Lip^(y)Lip(r/)||p||^. 
Proof. Let X,Y be the ODE flows associated to p, a. Let 70 G '^{p,a). Define 

Xt^Yf. {x,y) ^ {X{x),Y{y)). 
Then 74 = {Xf x lt)jj7o G 'B{pt,crf). Let us introduce 

Q{t)=l \x-y\d-ft{x,y)= I \Xt{x) -Yt{y)\d-iQ{x,y) . 

Then Q is a Lipschitz function and 

Q'{t)< [ \V{Xt{x),rt*ri{Xt{x)))-V{Yt{y),st*v{yt{x)))\djo{x,y). 

By triangular inequality, we obtain 

Q'{t) <(Lip,(y) + Lip,(y)Lip(n*r/))Q(t)+Lip,(y) / \{rt - st) * viYt{y))\ djo{x,y) • 
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Note besides that, thanks to Proposition 14.81 we have 

{ri-sl)*v{z)= [ v{z-0{drl{C)-dsi{0)<Up{v)W,(rlsi 



Integrating, we get 

Q{t) < Q{0)e^' + tC'e^'supWi{rt,st). 

T 

where C = Up^{V) + Lip^(y)||rt||_^Lip(r/), C = Lip^(y)Lip(7?)||p|l_;^^. 

We conclude taking 70 in an optimal way so that Q(0) = VVi(/00)Co) and using the in- 
equality 

m{pt,<Tt)<Q{t). 

□ 

The stability estimate allows us to apply Banach fixed point Theorem for T small enough. 
4.3 Measure solutions are Lagrangian solutions 

Proof of Proposition \4-6[ Let p be a measure solution of ([24|) . Let b = V{x,p * rj) and 



denote a the Lagrangian solution associated to dtcr + div(c75) = with a{0) = p. 

Then 6 = p — a is a measure solution of dtS + div{Sb) = 0, with S{0) = 0. That is to say, 
for any cp G - 00, T] x 



/ / (dt^ + b\t,x)-Vip)d5tdt = 0. 
Jo Jr^ ^ ' 



Let V G ^cQ - 'X),T\ X M^,M). We can find ^ G 'T^^Q - oo,r] x 1R^,R) so that = 
dt^p + b\t, X) ■ Vip. Hence, for any ^ € "^^G - oo,T] x R^, M), we have ijj d5t dt = 0, 

which implies 5 = a.e. so p = a a.e.. Consequently, we have U{t,x) = V'^{x,a * r/*), and p 
is a Lagrangian solution of ()24p . □ 



5 Conclusion 

In the Kruzkov framework, we are able to prove existence and uniqueness of weak entropy 
solution for the equation dtp + d\Y{pV{x, p, p * rj)) = 0. Furthermore, we can prove uniform 
bound in L°° \i V = v{p)W {x, p * rj), with v{l =)0. However, the required hypotheses are 
very strong: we need indeed V G n W^'^ n W^'°^. 

In the optimal transport theory framework, we can treat only equations such that dtp + 
d\Y{pV{x, p * rf)) = 0. For this equation, we have only L°° bound that are exponentially 
growing in time. The hypotheses are nevertheless weaker since we only ask V G Lip n L°°, 
but we are no longer able to prove the Gateaux-differentiability of the semi-group. 
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